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The SU(2) and SU(3) chiral phase transitions within Chiral Perturbation Theory 
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The SU(2) and SU(3) chiral phase transitions in a gas made of pions, kaons and etas are studied 
within the framework of Chiral Perturbation Theory. We describe the temperature dependence of 
the quark condensates by using the meson meson scattering phase shifts in a second order virial 
expansion. In particular, the SU(3) formalism yields a somewhat lower melting temperature for the 
non-strange condensates than within SU(2), and also predicts that the strange condensate melting 
is slower than that of the non-strange, due to the different strange and non-strange quark masses. 
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There is a growing interest in studying the phase dia- 
gram of QCD, namely, the transition from a hadron gas 
to a quark gluon plasma, the chiral phase transitions, 
their dependence on the number of flavors, on quark 
masses, possible color superconducting states, etc. 

The virial expansion is a simple and successful ap- 
proach to describe many thermodynamic features of di- 
lute gases made of interacting pions B and other hadrons 
0). Remarkably, for most thermodynamic properties it 
is enough to know the low energy scattering phase shifts 
of the particles, which could be taken from experiment, 
avoiding any model dependence. However, if one is in- 
terested in the chiral symmetry restoration, whose order 
parameters are the quark condensates, denned as deriva- 
tives of the pressure with respect to the quark masses, 
one needs a theoretical description of the mass depen- 
dence of the scattering amplitudes. 

With that aim we turn to Chiral Perturbation Theory 
(ChPT) H EL ||] , which provides a systematic description 
of low energy hadronic interactions. After identifying 
the pions, kaons and the eta as the Goldstone Bosons of 
the QCD spontaneous chiral symmetry breaking (pseudo 
Goldstone bosons, due to the small light-quark masses), 
the ChPT Lagrangian is built as the most general deriva- 
tive and mass expansion, over AttF ~ 1.2 GeV, (the 
symmetry breaking scale) respecting the symmetry con- 
straints. At one loop any calculation can be renormalized 
in terms of just a set of parameters, Lk(n), H^(^l), (/z be- 
ing the renormalization scale) which can be determined 
from a few experiments and used for further predictions 
at low temperatures. 

In particular, within SU(2) ChPT the leading coeffi- 
cients of the low temperature expansion for the pressure 
and the quark condensate (qq) = (uu + dd) in a hadron 
gas whose only interacting hadrons were the pions and 
the other hadrons were free, have been calculated j(| ob- 
taining also an estimate, T ~ 190 MeV, of the phase tran- 
sition critical temperature. It was also shown that the 
perturbative calculation of the pressure was analogous 



to the second order virial expansion when the interact- 
ing part of the second virial coefficient is obtained from 
the one loop 7T7t ChPT scattering lengths. Other works 
have studied the applicability of the virial expansion with 
a chemical potential 7], or the critical temperature in 
generalized scenarios [3 . 

The SU(2) approach is limited by the absence of other 
interacting particles, like kaons, whose densities at T ~ 
150 MeV are significant ||. Let us recall that the chiral 
phase transition can be different for the SU(2) and SU(3) 
cases, and that several QCD inequalities and lattice re- 
sults suggest a chiral condensate temperature suppres- 
sion with an increasing number of light flavors || , that 
is, entropy and disorder. This effect would lower the chi- 
ral critical temperature down by roug hly 20 MeV |l(J 
(in the chiral limit). Note that none of these results have 
been obtained from the hadronic phase. 

In addition, within SU(2) it is not possible to study 
the (ss) condensate. Furthermore, in the chiral phase 
diagram, quark masses play the same role as magnetic 
fields in ferromagnets: Intuitively, we need a higher tem- 
perature to disorder a ferromagnet when there is a mag- 
netic field aligned along the direction of the magneti- 
zation. Analogously, it has been found that the SU(2) 
chiral condensate melts at a lower temperature in the 
massless limit Q. Hence, since the strange quark mass 
m s > m u ,irid, we expect a "ferromagnetic" effect and a 
slower temperature evolution for (ss) than for (qq). Us- 
ing a QCD-motivated effective lagrangian in the large N c 
limit, this effect has been studied in Jll| 

In this work we will use SU(3) ChPT together with the 
virial expansion to estimate the effect on (qq) of pions in- 
teracting with kaons and etas, and study the effect of the 
different quark masses. In particular, we will obtain the 
temperature dependence of (ss) in the hadronic phase. 

Thus, the second order relativistic virial expansion of 
the pressure for an inhomogeneous gas made of three 
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species: i = w, K, rj is [El |i"2| 



/3P = 2 B<(T)6+£ ( B ug + \T, B iMi (1) 

where /3 = 1/T and ^ = exp(— flrrii). Expanding up to 
the second order in £j means that we consider only binary 
interactions. For a free boson gas, = for i ^ j, 



whereas: 



^ Jo 



(2) 



B <j>) = _£L / dpp 2 e- 2 ^V^f-™.), ( 3 ) 
4tt 2 7 

and the degeneracy is <?i = 3, 4, 1 for ir, K, r\, respectively. 
The interactions appear through || |l2|| : 



,(mt) _ Si 

y ~ 2 



/ dEE 2 K 1 (E/T)A^(E), (4) 



where is the first modified Bessel function and: 



A«= £(2/ + l)(2J+l)S? tJtS (E), 



(5) 



(5/* j s being the ij — * ij phase shifts (chosen so that 5 = 
at threshold) of the elastic scattering of a state ij with 
quantum numbers /, J, S (J being the total angular mo- 
mentum and S the strangeness). 

The virial expansion seems to break at T ~ 200 — 
250 MeV. Fortunately, the physics we are interested in 
occurs already at T < 250 MeV, and 4v > £k — 6yi so 
that it is enough to consider ij = irir, irK and irn in the 
second virial coefficients. Thus we are considering up to 
£ir suppressed corrections to the w, K, r\ dilute free gas. 

Let us recall that the pressure is nothing but the T 
dependent part of the free energy density z, i.e., P = 
eo — z, eo being the vacuum energy density. Hence, H: 



dz 



{q a q 



dm„ 



= (0|? Q ? Q |0) - 



dP 

dm 0n 



(6) 



where q a — u,d,s. Let us emphasize that, in contrast 
with most thermodynamic quantities, for the chiral con- 
densate it is not enough to know the 5(E), but we also 
need their dependence with the quark masses as well as a 
value for the vacuum expectation value. All that will be 
obtained from ChPT. For that purpose, we first rewrite 
the condensate, eq.(^), in terms of meson masses: 



{q a q a 



x q a \0) U + 



<?£><? a 



dP 



(T) 



2miF 2 drrii . 

where, as before, i = tt, K, rj. Note that we have defined: 



(8) 





Refs.[4 


, 22] 


Ref.[ 


21] 


IAM | 


L7 1 


L\(M P ) 
U 2 (M P ) 

L 3 
17AM \ 
Ll(M p ) 
L r 6 (M P ) 

L 7 
L%{M P ) 


0.4 ±0.3 
1.35 ±0.3 
-3.5 ±1.1 
-0.3 ±0.5 

1.4 ±0.5 
-0.2 ±0.3 
-0.4 ±0.2 

0.9 ±0.3 


0.46 
1.49 
-3.18 


1.46 



-0.49 
1.00 


0.561 ± 0.008 (±0.10) 
1.21 ±0.001 (±0.10) 

-2.79 ±0.02 (±0.12) 

-0.36 ±0.02 (±0.17) 

1.4 ±0.02 (±0.5) 
0.07 ±0.03 (±0.08) 

-0.44 ±0.003 (±0.15) 
0.78 ±0.02 (±0.18) 



TABLE I: Different sets of chiral parameters xlO 3 . In the 
second column L\,L2,Ls are taken from J22| and the rest 
from Q (L\ and L\ are estimated from the Zweig rule). The 
third column comes from an 0(p 4 ) analysis of Kn decays |5l]] 
(La and L% are set to zero). The last column is the IAM fit 
Wn to the whole meson-meson scattering up to 1.2 GeV. 



In the isospin limit, when m u — ma, the u and d conden- 
sates are equal, and we define (0|fg|0) = (0\uu + dd\0). 
It is tedious but straightforward to obtain expressions 
for the c parameters above. They can be found in the 
appendix. Both the (0|gq|0) and (0|ss|0) one loop cal- 
culations within SU(3) ChPT were given in ||, together 
with the one loop dependence of the meson masses on 
the quark masses needed to obtain dm 2 /dm qa . The only 
relevant comment is that the c coefficients depend on the 
chiral parameters Lk, for k = 4. ..8, and H2. 

In Table I, we show several Lk determinations from me- 
son data. It makes little difference to use other determi- 
nations. For H 2 we will use H%(M p ) = (-3.4± 1.1)10~ 3 , 
obtained as explained in ]l3| but using a more recent 
estimation of (0|ss|0)/(0|gg|0) =0.75 ±0.12 @. For in- 
stance, the c parameters obtained when using the Lk in 
the first column of Table I are: 



hO.4 



- qq — Q+°- 2 r qq — 5" 

-7r — u - J -0Ai c K — v -°-0.7i 



005+ 029 r Ss 

u.uuo_ 037 , C K 



1.3 



-0.4 
-0.8' 



0.4 
= f .5; 



+0.5 

-0.7' 

-0.9 
-1.6- 



Note that our c qq is in a good agreement with the SU(2) 
estimates: 0.85 and 0.90 ± 0.05 §. Had we used the 
estimate H%(M p ) ~ 2LZ(M p ) ~ 1.8 x 10~ 3 from scalar 
resonance saturation [fy| , all our following results would 
remain unchanged for (qq), but the strange condensate 
would not even seem to melt within this approximation. 
Finally, we need the meson-meson amplitudes to one-loop 
in SU(3) ChPT, which have been given in fT|, 0. These 
expressions have to be projected in partial waves of defi- 
nite isospin I and angular momentum J, whose complex 
phases are the 5ij(E) in eq.(|5|). ChPT is known to pro- 
vide a good low energy description (E < 500 MeV) of 
the meson-meson amplitudes |l6|. Within SU(2) ChPT 
it has been shown that using just the amplitudes at 
threshold (scattering lengths) in the virial expansion is 
equivalent to expanding the partition function to third 
order in T/F or T/m v . This approach yields a fairly 
good representation of the pion gas thermodynamics at 
low temperatures, T < 150 MeV M. 
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The virial coefficients and dP/drrii have been calcu- 
lated numerically. For the figures, we have represented 
the chiral condensate over its vacuum expectation value, 
(QaQa) I (0\<ja<la\ty j so that all of them are normalized to 
1 at T = 0. In reality, since there is always an small ex- 
plicit chiral symmetry breaking due to the quark masses, 
the condensate docs not vanish completely, except in the 
T — ► oo limit (following with our analogy, our ferromag- 
net above the Curie point becomes paramagnetic in the 
presence of a magnetic field). Of course, with the second 
order virial approach we cannot generate such an analytic 
behavior, and our curves become negative above some T. 
Since the non-strange quark masses are very tiny com- 
pared with the (qq) (about a 3%), it is a good approx- 
imation to let the condensate vanish and talk about a 
melting temperature. This is not the case for the strange 
condensate, which has a larger explicit breaking due to 
the much larger mass of this quark. 

In Figure 1 we show the results of using the one- 
loop ChPT expansion for the phase shifts. The thin 
continuous line represents the interacting pion gas in 
SU(2), but updated to the SU(2) parameters obtained 
from those in the first column of Table I (see |J for 
the conversion), whereas the thin-dashed line is the cor- 
responding one-loop SU(3) result also for a pion gas. 
We have checked that the tiny difference between them 
comes only from the 0(p 4 ) phase shift contribution, ( it 
is the same for SU(2) and SU(3) only up to sIM?^ or 
s /My terms from expanding kaon or eta loops Q ) . This 
amounts to a 4 MeV decrease of the pion gas extrapo- 
lated melting temperature: T^ q) = 231 MeV for SU(2) 
and T$ q) = 227 MeV for SU(3). Next, the thin-dotted 
line is obtained by adding free kaons and etas to the pion 
gas, which lowers the extrapolated melting temperature 
down to 221 MeV. 

Our new SU(3) result for a,ir,K,r] gas, is the thick con- 
tinuous line, where we see that the additional decrease 
due to irK,irr] interactions amounts to 10-11 MeV. It 
may seem surprising that the irK , tttj interaction effect 
may be comparable or larger than that of free kaons 
and etas, since they are thermally suppressed. Note, 
however, that exp(— m„/T) just amounts to a factor 
of 6,4,2.5 suppression at T= 80,100,150 MeV, respec- 
tively. In contrast, the irK and irrj interactions depend 
strongly on m v , much more sensitive to to than mx or 
nifj, which are the only dependence of the free K and r\ 
terms on to. In particular dm\jdm ~ 2dm 2 K /dm and 
dm^l dm ~ 6 dmx / 'dm, and this "temperature indepen- 
dent enhancement" competes with the thermal suppres- 
sion, making the irK and irrj interaction effect compara- 
ble to the free one around T — 80— 100 MeV and larger if 
T > 140 MeV. All together, the T r (qq) decrease in SU(3) 
is roughly 20 MeV, in agreement with the chiral limit lat- 
tice results T c = 173 MeV for SU(2) and T c = 154 MeV, 
for SU(3) 0. 



Furthermore, it can be noticed that (ss) melts much 
slower than (qq), since m s rh. Indeed, there is still 
70% left of the (ss) condensate at the (qq) melting point. 
Let us recall that both effects are already sizable at low 
temperatures T ~ 100 MeV, even lower when looking at 
the slower (ss) evolution, which already separates from 
the (qq) evolution at temperatures as low as 70 MeV. 
Since we can still trust pure ChPT up to, roughly, T ~ 
150 MeV §, we conclude that both effects are model 
independent predictions. 

However, just for illustration and to ease the compari- 
son between curves and with previous works, we have ex- 
trapolated the condensates to higher temperatures. We 
have also made a Montecarlo Gaussian sampling of the 
chiral parameters within their error bands. We find that 
the melting temperatures would be: 



T {q q )su(3) 

m 

T {qq)SU(2) 
m 

T (qq)SU(2) 



- 211I^ y MeV, 
: 231+?gMeV 

T (qq)SU(3) 



21+£ 4 MeV 



Note that the two melting temperatures are strongly cor- 
related. 

Concerning the (ss) extrapolation, its "melting" tem- 



perature T, 



(f S ) 



291±3lMeV, is just for illustration, 



since, apart from lying beyond the reliability region of 
the approach, as stated above, actually it should not melt 
completely, except in the T — > oo limit. In particular, we 
would find rif s> - T r {qq) = 80±H MeV. Nevertheless, we 
remark once more that the slower (ss) thermal evolution 
is clearly seen at temperatures where the whole approach 
can be trusted. As a matter of fact, from the figure, we 
see that the strange condensate does not start to melt 
sizably up to T ~ 150 MeV. 
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FIG. 1: Condensate evolution with temperature using one- 
loop ChPT amplitudes. Although they should not really van- 
ish they are extrapolated down to zero only for reference. 



Furthermore, we have also estimated the effect of 
other, more massive, hadrons, which in the SU(2) case 

H 0] also included the kaons and etas, and reduced Tm 9 ^ 
by approximately 10-20 MeV. Since they are heavier than 
to,;, their density is very low, and their main contribution 
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to the pressure comes from the first virial coefficient, i.e. 
the free gas. The only uncertainty is on dM^j dm qa , con- 
servatively estimated to lie within the number of valence 
quarks N qa and 2N qa . Thus, we have found that their 

contribution decreases Tm^ by 7-12 MeV, respectively. 
The heavy hadron contribution is the same for the SU(2) 
and SU(3) cases, so that the difference between their re- 
spective melting temperatures remains the same. Finally 
if we were to extrapolate the strange condensate evolu- 

(ss) 

tion we would also find a decrease of Tm by 15-23 MeV 
due to these heavier hadrons. 

As an attempt to estimate the uncertainty due to the 
truncation of the virial expansion, we have added to 
our calculation the third order coefficient for free pions, 
which can be calculated exactly. If we do so, the melting 
temperature decreases by one MeV. It seems that when 
higher orders corrections start becoming important, the 
curve is already so steep that the melting point changes 
very slightly. We cannot calculate the size of the three 
pion interaction term, but we expect it to be of the same 
order, and consistently with what we did for the heavy 
hadron uncertainty, we also allow for twice as much un- 
certainty, i.e., another 2 MeV. The important point, how- 
ever, is that the error is dominated by the uncertainty in 
the chiral parameters, which is an order of magnitude 
larger. Altogether, we find: 



Tjf> ~ 201+^ MeV 



(9) 



Finally, and in order to estimate the effects of higher 
energies we can extend ChPT up to E ~ 1.2 GeV by 
means of unitarization models 17, [j^, [lj|. These tech- 
niques resum the ChPT series respecting unitarity but 
also the low energy expansion, including the mass de- 
pendency. In particular, it has been shown that the cou- 
pled channel Inverse Amplitude Method (IAM) provides 
a remarkable and accurate description of the complete 
meson-meson interactions below 1.2 GeV, generating dy- 
namically six resonances from the one-loop ChPT expan- 
sion and a set of fitted Lk compatible with other previous 
determinations. This approach can be extended system- 
atically to higher orders (for SU(2) case, see po[). 

Hence, in Figure 2 we show the results of using the 
one-loop coupled channel IAM fitted phase shifts. Its 
corresponding L, parameters are given in the last col- 
umn of Table I, with two errors, the first, very small, is 
purely statistical, and the second covers the uncertainty 
in the parameters depending on what systematic error is 
assumed for the experimental data. Let us remark that, 
although the Lfc are highly correlated, the second, larger 
error, ignores completely these correlations, and should 
be considered as a very conservative range. The con- 
tinuous line corresponds to the central values, and the 
dark shaded areas cover the one standard deviation un- 
certainty due to the small errors in the parameters. These 
areas have been obtained from a Montecarlo Gaussian 



sampling. The conservative ranges are covered by the 
light gray areas. We now find: 



T {qq)SU(3) 
m 

T (qq)SU(2) 
m 

T {qq)SU(2) 



■■ 204^ MeV, 

: 235±? Ci ) MeV, 
T msu(3) = 31.50+1-20 ( 9 )MeV 



where the errors in parenthesis are the conservative er- 
rors which should be interpreted as uncertainty ranges 
better than as standard deviations. Note the excellent 
agreement with standard ChPT. The magnitude of the 
different contributions is roughly the same, although the 
ttK and irr) interactions this time lower Tm^ by 17 MeV, 
and the free kaons and etas by roughly 10 MeV. Let us 
remark that the strange condensate again melts much 
slower than the non-strange condensate, and that this 
effect is clearly sizable at low temperatures, where both 
ChPT and the virial expansion can be trusted. With 
all the caveats presented above, and just for illustration 



we also extrapolate (ss) to zero and we find: T, 



{ss) 



304^25 di) MeV, in particular, the strange condensate 



melting is retarded by T, 



(Ss) 



T^ q) = 100 



+36 (120 
-29 



OMeV. 



0.8 
0.6 
0.4 
0.2 



Temperature evolution 
of chiral condensates 

ChPT+IAM 




50 100 150 200 

Temperature (MeV) 



300 



FIG. 2: Temperature evolution of chiral condensates using 
the unitarized ChPT amplitudes. The shaded areas cover the 
uncertainties in different sets of chiral parameters. Although 
they should not really vanish they are extrapolated down to 
zero only to ease their comparison and for reference. 

Again other massive states apart from kaons and etas 
would lower T^ q) by 6-10 MeV, and T^ s) by 22-32 MeV. 
Also very similar is the effect of including the third or- 
der virial coefficient for free pions. As it happened in 
the standard ChPT case, it decreases the melting tem- 
perature by one MeV, and we add another 2 MeV of 
uncertainty for the unknown interacting part, that is 



195+J (V) MeV 



(10) 



Again, the errors correspond to the uncertainty in the 
chiral parameters plus that on the heavier states, plus 
that from the truncation of the virial expansion and 
those in parenthesis should be interpreted as conserva- 
tive ranges. 
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These results show that higher energy effects do not af- 
fect very much our conclusions and that the results of the 
integrals in eq.(Q), which extend to infinity, seem robust. 

We have studied the SU(2) and SU(3) temperature 
evolution of the chiral condensates in the hadronic phase. 
The thermodynamics of the meson gas has been obtained 
from the virial expansion and Chiral Perturbation The- 
ory. Our results clearly show a significant decrease, about 
20 MeV, of the non-strange condensate melting tempera- 
ture, from the SU(2) to the SU(3) case, similar to lattice 
results. Of these, about 6 MeV had already been ex- 
plained with free kaons and etas, but the rest are mainly 
due to ttK and nrj interactions. In addition, our results 
show an slower temperature evolution of the strange con- 
densate, about 90% of its zero value still remaining at 150 
MeV, due to the different quark masses. Both effects are 
clearly visible already at low temperatures. However, we 
have checked that their size is completely similar when 
using unitarized models. 

These techniques should be easily extended to Heavy 
Baryon Chiral Perturbation Theory, in order to study 
the condensates with non-zero baryon density. In gen- 
eral the whole approach could be used with any effective 
Lagrangian formalism, and in particular to study other 
QCD phase transitions like those of the color supercon- 
ducting phases. 

Appendix.- We provide here the expressions for the 
SU(3) a coefficients: 



vk. 



1 + 4/^ + 2fj, K + m 7r (z/ Tr — ) 



16m 



F 2 



K T 



4m 



F 2 



?(8L 6 + 6L S - 6L 4 - 4L 5 - H 3 ) 



„99 _ 
-K — 



1 

2 

2m 



3^7T 

2 



2 



Fo 2 



L 4 + L 8 + -^ 



Hj 
2 



8m 



f(4L 6 + 2L 8 -3L 4 -L 5 ) 



- jl + 6fiK - 3m^7r + y-( m ^ _ 4m 2 ) 



16m 



^0 

I6-L4 128 

ci2~ ( m K + m l/ 2 ) ~ TT^( m K ~ m D( 3L 7 + L *) 

^0 Mo 

-^L(2L 8 + H 2 )} 



F, 







16m 



^(2L 6 - U) - -v v mt 



4m 2 



{1 + 2 M „ + in K + — f [2(L 8 - L 4 ) + H 2 ] 



8 16m 2 

+ -Vr,m 2 K + — 3^(2L 6 + L s - 2L 4 - L 5 - H 2 /2)} 

9 ^0 



4 s = -{1 + 6/iK + m*{3v K /2 - 4^/3) 



-ml (v K j2 + j/,/3) 



16m 



K 



F 2 



4L 7 



L 4 - 



Lb 
3 



8m| 
F 2 



5 H 2 
—Li — -Lg + — 8X7 



m" ( 3i 4 
£ [ 4L 8 - 2L b - -A + 3L« 



where following the standard notation in rcf.|^| 
I* = mflog(m 2 / / i 2 )/(32 7 r 2 J F 2 ), v = (log^ 2 /^ 2 ) + 
1)/(327t 2 Fq) and Fo is the pion decay constant in the 
SU(3) chiral limit. In this work we have set the renor- 
malization scale to \x = M p 
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